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Abstract 

The joint eigenvalue distributions of random-matrix ensembles are de- 
rived by applying the principle maximum entropy to the Renyi, Abe and 
Kaniadakis entropies. While the Renyi entropy produces essentially the 
same matrix-element distributions as the previously obtained expression 
by using the Tsallis entropy, and the Abe entropy does not lead to a 
closed form expression, the Kaniadakis entropy leads to a new general- 
ized form of the Wigner surmise that describes a transition of the spacing 
distribution from chaos to order. This expression is compared with the 
corresponding expression obtained by assuming Tsallis' entropy as well as 
the results of a previous numerical experiment. 

1 Introduction 

In 1957 Jaynes [I] proposed a rule, based on information theory, to provide 
a constructive criterion for setting up probability distributions on the basis of 
partial knowledge. This leads to a type of statistical inference which is called the 
maximum-entropy principle (MaxEnt). It is the least biased estimate possible on 
the given information. Jaynes showed in particular how his rule, when applied to 
statistical mechanics, leads to the usual Gibbs' canonical distribution. The core 
of the MaxEnt method resides in interpreting entropy, through the Shannon 
axioms, as a measure of the "amount of uncertainty" or of the "amount of 
information that is missing" in a probability distribution. This was an important 
step forward because it extended the applicability of the notion of entropy far 
beyond its original roots in thermodynamics. In this paper we consider the 
application of MaxEnt to the random-matrix theory (RMT), which is often 
used to describe quantum systems whose classical counterparts have chaotic 
dynamics [21 El- This is the statistical theory of random matrices H whose 
entries fluctuate as independent Gaussian random numbers. Dyson 4 showed 
that there are three generic ensembles of random matrices, defined in terms of 
the symmetry properties of the Hamiltonian. The most popular of these is the 
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Gaussian orthogonal ensemble (GOE) which successfully describes many time- 
reversal-invariant quantum systems whose classical counterparts have chaotic 
dynamics. Balian 5; derived the weight functions P{H) for the random-matrix 
ensembles from MaxEnt postulating the existence of a second moment of the 
Hamiltonian. He applied the conventional Shannon definition for the entropy 
to ensembles of random matrices as 
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J dHP{H)h\P{H) (1) 



and maximized it under the constraints of normalization of P(H) and fixed 
mean value of Tr(iJ T if), where the superscript T denotes the transpose. The 
latter constraint ensures basis independence, which is a property of the trace of 
a matrix. Then, the distribution P{H) is determined from the extremum of the 
functional 

Fsch = S Sch -iJdH P(H) - V JdH P(H)Ty (H t H) , (2) 

where £ and rj are Lagrange multipliers. Its maximum is obtained equating its 
functional derivative to zero. He obtained 

P{H) cc exp [-r)Tr (H T H)] (3) 

which is a Gaussian distribution with inverse variance l/2?7. Most of the in- 
teresting results are obtained for N x N matrices in the limit of N — > oo. A 
complete discussion of the level correlations for large matrices is a difficult task. 
Analytical results have long ago been obtained for the case of N — 2 The 
distribution of nearest-neighbor-spacings (NNS) of levels for GOE of 2 dimen- 
sional matrices, 

P(s) = (vr/2) s exp (-7rs 2 /4) , (4) 

is known as Wigner's surmise. Here s is scaled so that the mean spacing equal 
one. Analogous expression are obtained for the Gaussian unitary and symplectic 
ensembles. To demonstrate the accuracy of the Wigner surmise, we expand this 
distribution in powers of s to obtain 

p( s ) = 32 2 (l - i s 2 + . . . \ 3 242s 2 _ 4 128s 4 + . (5) 

while the power-series expansion of the corresponding exact distribution 
yields 

7T 2 27T 4 

fcxact(s) - ^rs 2 - —s i + ■■■ = 3.290s 2 - 4.329s 4 + • • • . (6) 
3 45 

The Wigner surmise has been successfully applied to the NNS distributions for 
numerous chaotic systems. 

For most systems, the phase space is partitioned into regular and chaotic 
domains. These systems are known as mixed systems. Attempts to generalize 
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RMT to describe such mixed systems are numerous; for a review please see [3]. 
MaxEnt can contribute to the generalization of RMT by either introducing ad- 
ditional constraints or modifying the entropy. An example of the first approach 
is the work of Hussein and Pato [7], who use MaxEnt to construct "deformed" 
random-matrix ensembles by imposing different constraints for the diagonal and 
off-diagonal elements. 

The second approach to include mixed systems in RMT is to apply MaxEnt 
to entropic measures other than Schannon's. Recently, the Tsallis statistical 
mechanics [HI E] has been applied to include systems with mixed regular-chaotic 
dynamics in a nonextensive generalization of RMT |1UI 111! IT21 1131 IT^lll5| . Dif- 
ferent entropic forms have been introduced, like Tsallis', by appropriately de- 
forming the logarithm in the expression for the Shannon entropy |16j. Among 
the many generalizations, one can find the entropies considered by Renyi |17|. 
by Abe ^BEI|> by Landsberg and Vedral [23! , an d by Kaniadakis (re-entropy) 
PTl [221 EI] ■ The present paper is a natural continuation of the previous work 
on nonextensive RMT, by considering other forms of generalized entropy. 



2 Nonextensive RMT 

In this section, we review the main results of non-extensive RMT based on 
Tsallis entropy, and show that replacing it with other generalized entropies 
does not always produce new results. The matrix-element distribution P{H) is 
obtained by extremizing Tsallis' g-entropy 



S Ts [q, P{q, H)} = 1 - / dH [P Ts (q, H)] q / (g - 1) 



(7) 



rather than Shannon's, but again subject to the same constraints of normaliza- 
tion and existence of the expectation value of Tr (i/ T i/). This can be done by 
replacing 5g c h m Eq. (2) by Sts and then equating the functional derivative in 
the resulting expression to zero. The extremization then yields 



PTs(q,H) 



Z- 1 [l + (q-l) Vq Tr(H T H)]- 



(8) 



where Z q and r\ q are expressed in terms of the Lagrange multipliers. The case of 
2-dimcnsional ensembles belonging to the three symmetry classes is considered 
in Ref. For ensembles with orthogonal symmetry, the resulting NNS 

distribution is given by 



Pts(s) 
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-l/(9-l) + l/2 



(9) 



for q > 1, where 



r[(i/( g -i)-2] 
r [i/(g-i) -3/2] 



(10) 
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and r[x] is Euler's gamma function (23- The extension of these results to the 
case on ensembles with arbitrarily large N amounts to replacing q in Eq. (9) 
by another jV-dependent parameter |15) . 

It is well known that the Shannon entropy of the couple of two independent 
random variables is the sum of their respective Shannon entropies. We have 
firstly attempted to extremize Renyi's entropy for the matrix-element probabil- 
ity density function P Rn (H), which is defined by 

S Rn [a, P Rn (H)} = ^— In (^J dH [P^{H)] a \ , (11) 

where the index a varies between 1 and 2. The obtained distribution NNS is very 
similar to that in Eq. (9), which is obtained by the maximization of Tsallis' 
entropy, if we replace 1 — q by a — 1 . In other words, mere fittings of observed 
data by the NNS distributions in Eq. (9) do not tell us anything about which of 
Tsallis' or Renyi's entropy is the underlying physical quantity. This result agrees 
with the known fact |25| that the maximum entropy principle in combination 
with Renyi's entropy reproduces the equilibrium probability distributions of 
Tsallis' non-extensive thermostatistics. In fact, the Renyi entropy is, in all 
cases, a monotonically increasing function of the Tsallis entropy 



ln[l + (l-g)P Ts ] 



1-9 



(12) 



Therefore the probability density function that extremizes one of these entropies 
automatically extremizes the other. 

Abe's entropy |^] is another interesting example of non-extensive entropies. 
It applies a g-calculus, which is invariant under the exchange q — > q^ 1 . This 
is a symmetry which plays a central role in the physical context of quantum 
groups 27 . This entropy has been introduced in |181 ITU] and has been applied 
there to generalized statistical-mechanical study of g-deformed oscillators. Abe's 
entropy, generated for the distribution function of a random-matrix ensemble, 
is defined by 

S Ahc [q, P(q, H)] = - j dH ([P A b.(«, H)] 9 - [P Abe (g, H)]~ q ) j (q ~ q' 1 ), 

(13) 

and also reduces to Schannon's entropy in the q — > 1 limit. Due to the q — > q 1 
symmetry the range of q may be restricted to < q < 1. Now, replacing Ssch 
in Eq. (2) by »!?Abc an d then equating the functional derivative in the resulting 
expression to zero, one obtains after some reorganizations 



^be - ^AbcTr (H T H) = 0. (14) 



q-q- 



Unfortunately, it is not possible to solve this equation with respect to i"Abe 
analytically. The analogous equation is solved numerically in Ref. ^Sl- When 
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q = 1, one obtains the standard GOE distribution, whereas for < q < 1, the 
distribution comes to exhibit the power-law behavior similar to similar to the 
distributions obtained by extremizing the Tsallis and Renyi entropies. Indeed, 
in the asymptotic region where i"Abc *C 1, one neglects the term in which PAbe 
is raised to a positive power and obtains 

P Ahc (q, H) ~ [£ Abc + 77 Abc Tr (H T H)} . (15) 



3 Kaniadakis' entropy 

In this section, we consider a possible generalization of RMT based on an ex- 
tremization of Kaniadakis' K-entropy |21U22ll2*5| . This entropy shares the same 
symmetry group of the relativistic momentum transformation and has appli- 
cations in cosmic-ray and plasma physics. For the matrix-element probability 
distribution function, it reads 

5 K [«, PkK H)] = ~ fdH (-^- [Pk(k, H)] 1+k - ^— [P k (k, H)} 1 ^ 
Ik J \ 1 + k 1 — k 

(16) 

with k a parameter with value between and 1; the case of k = 0, corresponds 
to the Schannon entropy. Here, a is a real positive parameter. Kaniadakis 
has considered two choices of a, namely a = 1 and a = Z, where Z is the 
generalized partition function. We here adopt the second choice. The matrix- 
element distribution Pk{k,H) is obtained by extremizing the functional 

F k = S k -VkJ dH P k (k, H)Tt (H t H) , (17) 

where r] K is a Lagrange multiplier. One arrives to the following distribution 

P k (k, H) = | exp {K} [-r^Tr (H T H)] , (18) 

where 

Z = JdH exp {K} [-7j K Tr (H T H)] . (19) 
Here exp^ K j. [x] is the K-deformed exponential j21j which is defined by 

expj K j [x] — + k 2 x 2 + x^j — exp ^—arcsinh kx^J . (20) 

and has the properties exp.r } [x] = exp(cc), exp^j [0] = 1 and exp|-_ K j [x] — 
exp^ K i [x] , and obeys the scaling law exp,r K j. [Ace] = (exp^ K ^ A j. [x]) A .The asymp- 
totic behavior of the K-exponential is expr K i [x] ~ |2kcz;| ±1 /' k ' . 

* ' x — >±oo 

To calculate Z, we note that Tr (H T H) = Yn=i H l + 2 J2 t>3 H ?r We 
introduce the new coordinates r = {r±, ■ ■ ■ ,r n }, where n = N(N + l)/2, and 
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rf stand for the square of the diagonal elements or twice the square of the non- 
diagonal elements, respectively. With these new variables, Eq. (19) becomes 



z = 2 -N(N-l)/4 / d n rexp 



-—arcsinh Krj K r' 



-N(N-l)/4 



ra/2 



T (— — -) 



(i + W2)r(^ + f) 



(21) 



for k < 2/n. Here we use the result obtained by Kaniadakis in evaluating an 
analogous integral that appears in his treatment of Brownian particles (Eq. (71) 
ofRef. ED)- 

We now calculate the joint probability density for the eigenvalues of the 
Hamiltonian H . With H = U~ 1 EU, where U is the global unitary group, we in- 
troduce the elements of the diagonal matrix of eigenvalues E = diag(£a, ■ ■ • , En) 
of the eigenvalues and the independent elements of U as new variables. Then 
the volume element dH has the form 



dH = \A N (E)\dEd(j,(U), 



(22) 



where Ajy {E) — Yl n>m {E n — E m ) is the Vandermonde determinant and dfJ,(U) 
the invariant Haar measure of the unitary group [3 El- The probability density 
Pk(k, H) depends on H through Tr(fftf/J and is therefore invariant under 
arbitrary rotations in the matrix space. Integrating over the " angular variables" 
U yields the joint probability density of eigenvalues in the form 



P k (k; E lt E N ) =C Y[(E„ -£ m )exp {re} 



N 
i=l 



(23) 



where C is a normalization constant. 

In order to obtain a generalization of the Wigner surmise, we consider the 
case of two-dimensional random-matrix ensemble where N — 2 and n = 3. In 
this case, Eq. (23) reads 



Pk(k;£,s) = 
where e = (E\ 



2(1 + 3k/4) 



B 



, 2k 



3 cP 
4' 2l 



(K?7 K ) 3/2 sexp {K} 



"?7k 



2 s 



1 , 



(24) 

Y{a)Y{b)/Y{a + b) is 



E 2 )/2, s = [Ex-Ei], and B{a.b) 
the Beta function [21]. The NNS distribution is obtained by integrating (24) 
over e from — oo to oo. This can be done by changing the variable e into 
x = exp[— iarcsinh(/i7/ K s 2 /2)], integrating by parts, and then replacing the 
variable x by another new variable, y — exp(«;a;). The resulting integral can be 
solved by using the following identity (21] 



y-\y + P) v (y - nf' 1 dy = u^- x b (\-»-u,h) 



/3 N 



2Fi[-v,X-n-v;X-v;-^J, (25) 
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for \(3/u\ < 1 and < \i < A — v, where 2F\(v , /i ; A; x) is the hypergeomet- 
ric function. Thus, after straightforward calculations we can express the NNS 
distribution as 

P K (s) = -2 (l + ^ rfose^-V^arcsinhC^/a) Bij-j, 5 ) 

2^! f-J, - - 1 ; - + 1; _ e -2arc S i„h( KI)KS 2 /2)"j _ ( 26 ) 

The condition of unit mean spacing defines the quantity rj K as 

irk 3 / 2 (1 + |«) 



»?K = 



(27) 



We note that the function £? — 5,5) diverges at k = k c = 1/2, which serves 
as an upper limit for the range of variation of n. We also note that the mean 
square spacing diverges at k = =2/5. 

We now compare the spacing distributions Pk(s) and Pt s (s). Each of the 
the two distributions coincides with the Wigner distribution for the smallest 
values of the entropic indices k ~ and q = 1, but essentially differs from 
the Poisson distribution and practically behaves as a delta function S(s) at the 
maximum allowed values of k = 0.5 and q = 1.5. Therefore, neither Pk(s) nor 
Pts(s) can describe the initial stage of the transition from chaos to integrability. 
In fact, it is hard to believe that nearly integrable systems can be described 
by a base-independent random-matrix model like the one considered in the 
present paper. Both distributions increase linearly with s near the origin. The 
Kaniadakis distribution Pk(s) decays at large s as s 2 " 2 / K . The NNS in Eq. (9) 
obtained by the Tsallis statistics asymptotically behaves as Pr s (s) ~ s 2-2 /^- 1 ). 
Both distributions have similar asymptotic behavior when k = q — 1. However, 
distributions Pk(s) and Pr s (s) satisfying this relation have different behavior 
at intermediate values of s. The evolution of Pk(s) as k varies from the GOE 
value of to ~ 0.5 and Pts(s) as q varies from the GOE value of 1 to ~ 1.5 
is demonstrated in Fig.l. The figure shows two sets of NNS distributions, one 
set for each statistics, with entropic indices related by n = q — 1 so that the 
corresponding distribution have similar asymptotic behavior. While the peak 
of the Tsallis distribution steadily moves towards the origin as q increases, the 
peak of the Kaniadakis distribution remains almost in the same position as k 
increases from and only starts to move towards smaller s when n exceeds a 
value of ~ 0.4. We remind that the mean square spacing diverges in this latter 
domain of k. 



4 Comparison with numerical experiment 

This section offers a comparison of the generalized expressions of the Wigner sur- 
mise derived from the Tsallis and Kaniadakis entropies with the NNS obtained 
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in the numerical experiment done by Zyczkowski and Kus [2E|- This experi- 
ment imitates a one-parameter (denoted by 5) transition between an ensemble 
of diagonal matrices with independently and uniformly distributed elements, 
corresponding to 5 = 0, and a circular orthogonal ensemble for 5 = 1. These 
authors were able to achieve reliable statistics by constructing numerically 500 
matrices of size N = 100 for each considered value of the transition parameter 
5. Diagonalizing these matrices yielded NNS distributions consisting of 50 000 
spacings for each value of 5. Three of these distributions intermediate between 
the Wigner and Poisson distributions, corresponding to 6 = 0.9, 0.5 and 0.1, 
are reported in Ref. 28]. These distributions are compared in Fig. 1 with the 
spacing distributions in Eqs. (9) and (26), which have been obtained by using 
the Tsallis and Kaniadakis statistics, respectively. The figure suggests that both 
entropies yield reasonable description of the stochastic transition although the 
agreement with the numerical results is not complete. The slight difference in 
the behavior of the two non-extensive NNS distributions is a positive aspects 
since it is well known that the transition from order to chaos does not neces- 
sarily proceed through the same passage. It remains an open question to find 
which systems are better described by either entropy. 

5 Conclusion 

Tsallis' entropy has been considered by several authors as a starting point for 
constructing a generalization of RMT. However, there are several other gen- 
eralized entropies. This paper derives random-matrix ensembles by maximiz- 
ing the generalized entropies proposed by Renyi, Abe and Kaniadakis under 
the constraints on normalization of the distribution and the expectation value 
of Tr(77 T i7). The distribution functions obtained for the Renyi entropy are 
analogous to the corresponding results, previously obtained assuming Tsallis' 
entropy, which is not a surprise since these entropies are monotonic functions of 
each other. The use of Abe's entropy did not lead to closed-form expressions for 
the matrix-element distribution. On the other hand, we obtain a new analyti- 
cal formula for the distribution function in the matrix-element space when the 
entropy is given by the Kaniadakis measure. This leads to a new expression for 
the NNS distribution of the eigenvalues in the special case of two-dimensions 
random-matrix ensembles. This special case is the one that leads to the Wigner 
surmise when the entropy is given by the Schannon measure. The high accuracy 
of Wigner's distribution describing chaotic systems justifies the use of the two- 
dimensional ensembles for obtaining NNS distributions using other entropies. 
As in the case of Tsallis' entropy, the NNS distribution obtained here for Kani- 
adakis' entropy has an intermediate shape between the Wigner distribution and 
the Poissonian that describes generically integrable systems. However, neither 
of these distributions reaches the Poissonian form for any value of the entropic 
index. It is generally believed that the route from integrability to chaos is not 
unique. Therefore, we hope that the NNS distribution obtained in this paper as 
well as the one, previously obtained with the Tsallis statistics, are good candi- 
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dates for modelling systems with mixed regular-chaotic dynamics at least when 
they are not far from the state of chaos. We have tested these expressions by 
comparing their predictions with the NNS distributions obtained in a numerical 
experiment by Zyczkowski and Kus. 
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Figure caption 

FIG. 1. NNS distributions obtained by using the Tsallis (lower panel) and 
Kaniadakis (upper panel) statistics. The entropic indices used, being q = 1 
(GOE), 1.1, 1.2, 1.3, 1.4, 1.45, 1.47, 1.49 in the case of the Tsallis statistics and 
K = (GOE), 0.1, 0.2, 0.3, 0.4, 0.45, 0.47, 0.49 in the Kaniadakis', are chosen so 
that the corresponding curves in the two statistics have the same asymptotic 
behavior. 

FIG. 2. NNS distributions obtained in a numerical simulation [2Hj of the 
transition from the GOE statistics to the Poissonian (histograms) compared to 
NNS distributions obtained by using the Tsallis (dashed lines) and Kaniadakis 
(solid lines) statistics. 
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